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Abstract 

^..^ We introduce a CP violating coupling between the Ricci curvature and the heavy right handed 

^^ . neutrinos. This splits the Majorana masses of the left and the right handed heavy neutrinos. When 






o 



X 



the heavy neutrinos decay, their decay rates are different into left and right chirality light neutrinos. 
A time varying non-zero Ricci curvature can give rise to a net lepton asymmetry. The source of a 



^S^ • non-zero curvature in a radiation era is the quantum fluctations of primordial fields. 

^' PACS numbers: 98.80.-k, 98.80.Cq 
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1. INTRODUCTION 

The origin of the baryon number asymmetry in the Universe is, as well known, a still 
open problem of the particle physics and cosmology ll| . The successful prediction of the 
light element abundances by big bang nucleosynthesis [2, y] and the observations of Cosmic 
Microwave Background anisotropics combined with the large structure of the Universe [4, |5| 
show that the baryon to photon number of the universe, i.e. the parameter that characterize 
such a a asymmetry, is of the order 

r] = ^^^^<(9.2±0.5)10-^S (1) 

where hb {^3) is the baryon (antibaryon) number density, and s the entropy of the Universe 
(s = 27r^(7*sT'^/45, g^:s counts the total degrees of freedom for particles that contribute to the 
entropy of the Universe, and finally g^,s takes values very close to the total degrees of freedom 
of effective massless particles f?*). Explaining tib in terms of known parameters is one of 
the major problems in cosmology. In a (CPT invarinat) theory, the baryon asymmetry can 
be explained provided that the Sakharov' conditions are satisfied |6[: 1) there must exist 
processes that violate the baryon number; 2) the discrete symmetries C and CP must be 
violated; 3) departure from thermal equilibrium. However, the Sakharov's conditions may 
be relaxed in some circumstances p]. As shown in [sj, in fact, a dynamical violation of 
CPT (which implies a different spectrum of particles and antiparticles) may give rise to the 
baryon number asymmetry also in a regime of thermal equilibrium (see also (9|). 

In this paper raise the possibility that the baryon/lepton asymmetry in the Universe can 



be generated by a CP violating curvature coupling of the heavy Majorana neutrinos 10| 



£^ = V^/3RiJz^,ij, (2) 

where a non-zero Ricci curvature is generated in the radiation era by back-reaction effects 
of quantum fields. The dimensional /3 will be assumed to be of order Mp , where Mp = 
G^^'"^ is the Planck mass. Quantum effects, in particular at the very early time, cannot be 
ignored because they may modify the dynamics of the Universe evolution. General relativity, 
therefore, requires some modifications in order to account for them. In absence of a complete 
quantum field theory of gravitation, one works assuming a semiclassical theory of gravity. 



In this context, the Einstein field equations read [ll|, [l2| 

i?.. - \9.^R = ^ w + (nr^)) (3) 

where T^t is the stress energy-momentum tensor for the classical field, T^u represents the 
energy momentum tensor operator for quantum fields, and finally {T/hJ ) = {0\T^u |0) 
represents the regularized expectation value of T^u . During the radiation dominated era, 
the trace of the classical energy momentum tensor vanishes, T*^^'^ = 0. However, owing to the 
presence of the quantum corrections (0|T/jt; |0), it follows that the trace is nonvanishing. 
This anomaly comes from the infinite counterterms that must be add to the gravitational 
action to make the trace finite, and it is responsible, as we will see, of the generation of the 
gravitational baryon asymmetry. 

The lepton-asymmetry gets frozen-in when the GUT scale lepton-number violating pro- 



cesses decouple. Baryon asymmetry can then 
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3e generated from this lepton-asymmetry 



- Il5| by the electro-weak sphaleron processes 



16| . Sphaleron processes conserve {B — L) 



so a lepton asymmetry generated in the GUT era can be converted to baryon asymmetry of 



the same magnitude |l3 |. 



2. CP ODD GRAVITATIONAL COUPLING OF MAJORANA NEUTRINOS 

We consider the CP violating interaction between fermions and the Ricci curvature de- 
scribed by the dimension- five operator ([2]) [10|. This operator is invariant under Local 
Lorentz transformation and is even under C and odd under P and conserves CPT. In a 
non-zero background i?, there is an effective CPT violation for the fermions. The Dirac 
equation 

27^(9^7/^ -Mi)- i(3R 75^^ = , (4) 

leads to the dispersion relation 

S^V = (p' + M2 + /32i?2)^ _ [3{^^Yd^R)ilj . (5) 

For a spatially homogenous background, the energy levels of the left and right handed 
fermions split as 

EV± = (p' + M^ + (3^R^ ± l3R)i,± , (6) 



where iIj± = (1/2)(1±75)'?/' and over-dot represents time derivative. Now consider a Majorana 
neutrino in the chiral representation Nm = {NR,ia2N^)'^ = {Nr, Nr')'^ . R ^ imphes a 
spontaneous CPT violation. The energy levels of the right chirality neutrinos Nr = ip^ and 
the left chirality antineutrinos {i^lY = V"- ^ire 

E± = Vp2 + M2 + pm^ ± ^^ = . (7) 

The effective masses M± = E±{p = 0) of the left and the right chirality states are different, 

M±c^M+ ^4r ±^R (8) 

2M 2M ^ ^ 

This mass splitting between the heavy neutrinos and antineutrinos results in the decays. 

The heavy right-handed Majorana neutrino interactions with the light neutrinos and Higgs 

relevant for leptogenesis, are described by the lagrangian 

£ = -/i„^(0t W^Ilp) - ^N%MNr + h.c. , (9) 

where M is the right handed neutrino mass-matrix, l^^ = (^c^e")!^ is the left-handed 
lepton doublet (a denotes the generation), (p = (0+,0°)-^ is the Higgs doublet. In the 
scenario of leptogenesis introduced by Fukugita and Yanagida, lepton number violation is 
achieved by the decays Nr -^ (p + II and also Nr^ — )■ 0^ + II". The difference in the 
production rate of l^ compared to /£, which is necessary for leptogenesis, is achieved via 
the CP violation. In the standard scenario, n{NR) = n{N'^) as demanded by CPT, but 
T{Nr — !■ /l + 0) 7^ ^{^R ~^ ^L + 0^) due to the complex phases of the Yukawa coupling 
matrix hnB, and a net lepton number arises from the interference terms of the tree-level and 

nn 

one loop diagrams [M, 115| . 

In our leptogenesis scenario there is a difference between the heavy right and left chirality 
neutrinos at thermal equilibrium due to the CP violating gravitational interaction ([2]). Now 
even in the absence of CP violating phases in the couplings, there will be a net lepton 
asymmetry generated from the decay of unequal number of left and right chirality heavy 
neutrinos. The decay rate of a heavy neutrinos is 

r± = ^/i'M± (10) 

The lepton asymmetry will be given by 

_ r+ - r_ _ M+ - M„ _ /3i^ 
'^ r, +r_ M M2 ^ ^ 



The lepton asymmetry will freeze-in at the temperature where the GUT scale lepton number 
violating processes decouple. In the next section we show that a time varying Ricci curvature 
can be generated by quantum fluctuations of the flelds in the radiation era. 

3. TIME VARYING RICCI CURVATURE FROM QUANTUM FLUCTUATIONS 

The dynamical evolution of the gravitational background is assumed to be described by 
the FRW metric with curvature K = 0, i.e. ds"^ = dt^ — a^{t){dx'^ + dy'^ + dz^). The form of 
the scale factor a{t) is obtained from the semiclassical Einstein equations ([3]). 

The regularized components of the energy-momentum tensor we concern have the form 

HQ 

(Tir^) = A:i«i/,. + A:3(^)//,., (12) 



where 



p2 

^^^H^, = 2R,^,, - 2g^,UR + 2RR^, - —g^, , (13) 
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n = V^V^, and ; stands for covariant derivative. The coefficients fci^s are constants. No- 
tice that these coefficients come from the regularization process and their values strictly 
depend not only on number and types of flelds present in the Universe, but also on the 
method of regularization. Because the methods of regularization affect the the values of 
fci 3 and more important because of the uncertainty of what flelds were present in the very 
early Universe, they can be considered as free parameters |l2l . Il8| . A comment is in or- 
der. The tensor ^^^H^^, satisfles V^^^^Hj^ = and is derived by varying the local action: 

A r 
^^^Hnu = 2\/ —g- — / d^\/—gPc'. To cancel the inflnities in IT^^^'^^) one has to add inflnite 

ogi^v J 
counterterms in the Lagrangian density describing the gravitational flelds. One of these 

counterterms if of the form ^—gCR?^ and due to (the logarithmically divergent) constant 



C, the coefficients k\ can take any value, and can be fixed experimentally 19|. As regards 
^^^Hf^i^, it is covariantly conserved only for conformal flat spacetimes, and cannot be derived 
by means of the variation of a local action, as for ^^^H^^. The coefficient k^ takes the form 



h = o \No + — A^i/2 + 31A^i ) 

^ 14407r2 V 2 ^ J 



For a SU{5) model, for example, the number of quantum fields take the values A'^o = 34, 



A^- 



1/2 



45, and A^i = 24, so that h ~ 0.07 
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The explicit expression of the components of ^^^H^^ and ^^^H^j^^ are 
-noo 






(^)i7, 



18{2HH + H^ + 10HH^ 
6 I 2 



(14) 



00 



3H\ ^'^H, 



12HH + UHH^ + 7H' j gij 
H\4H + 3H^)g,j. 



00 



The process of regularization leads to the a trace anomaly 20 1 



(15) 



that for the FRW Universe reads 



(T(«*^)) = 36fci 



d^H 



+ IHH + 4if 2 + I2if if 2 ^ 



+ l2hH'^ [H + H 



(16) 

H = a/ a is the Hubble parameter. The energy-momentum tensor of classical fields is the 
usual of perfect fluid T^]^ =diag (p, — p, — p, — p), where p is the energy density, p the pressure. 
They are related by the relation p = wp, being w the adiabatic index. During the radiation 
dominated era, the equation of the state is p = p/3, i.e. w = 1/3. This implies that the trace 
of the energy-momentum tensor of classical fields vanishes T^'^''^ = p — 3p = 0. However, the 
anomaly trace gives a non vanishing contribution. To evaluate it, we recall that during the 
radiation dominated era the scale factor evolves as a{t) = (aot)^/^, whereas the cosmic time 
and the temperature are related by the relation -^ = ^^qq^* ^ ■ These expressions should 
be in principle modified by the backreaction effects induced by quantum fields, but as we 
shall see (and for our purpose) the evolution of the Universe can be described by standard 
cosmology. 

For the FRW, the modified Einstein field equations assume the form 



3W 
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p + 18ki{2HH + H^ + lOHH^) + 3hH^ 



Ml . 



3H^ + 2H 



»7r 
Ml 



p + QkA2 



d^H 



dt^ 



(17) 
(18) 



+12HH + UHH^ + 7H^) + k^H^{4:H + 3H 



from which it follows 




2H'^ + H = 


Stt 
Ml 



Qhi^T- + ^HH + AH'^ + 12HH'^) 
dt-^ 

+ 2hH\H + H^)\. (19) 

We are looking for solutions of the form 

H{t) = H^{t) + 5{t), (20) 

where 5{t) ^ 1 is a perturbation, and Hq = l/2t is the Hubble parameter for a Universe 
radiation dominated. During the radiation dominated era, the Ricci curvature vanishes, 
i? = 0, as well as its covariant and (cosmic) time derivatives (Di? = 0, V^V^R = 0, R = 0, 
R = 0). It then follows that ^^'^H^^{H = Hq) = (i.e. ^ + 7HoHo + AH^ + 12HoHl = 0) so 
that it only contains 5-terms and its derivatives. On the contrary, ^'^■'if^jjy 7^ when H = Hq 
owing to the (quadratic) Ricci-terms (in other words we have Hq{Hq + Hq) ~ t~^ + Fs, 
where Fs is a function of 6 and its derivatives). This implies that in the final expression of 
the baryon asymmetry rjL, only the /cs-terms appear. 

Inserting H given in ( 120|) into Eq. ( TT9|) . to leading order one obtains 

• 2 ^ 167rA;3 1 ^ _^. 

whose solution is 

Sit) ~ Al _ ^Ll (22) 

C is a constant of integration. As it can be seen, the Mp'^ suppresses considerably the effects 
of 5 on the dynamics of the Universe evolution, and these terms wash-out for large t. During 
the radiation era H ~ l/2t, it follows that the trace anomaly (1T6!) reads 



(XiQM)'^ 



3fc3 



Stt 



From R = -772" (T^'^*^)) we find that 



Ml 



Ml t^ ^ ^ 

the parameter characterizing the heavy neutrino asymmetry fITT]) assumes the form 



,.(r)^*32W^^-|53^.3M,10»-^-^. (24) 



where 13 ~ Mp^. The heavy neutrino asymmetry freezes at the value rfLiTD) where Td is 
the decouphng temperature when the lepton-number violating interactions {Nr ^ N'^) go 
out of equilibrium. In the following section we will calculate To and rjiiTj:,) in the context 
of a specific GUT model as an example. 

4. HEAVY NEUTRINO ASYMMETRY AT THERMAL EQUILIBRIUM AND 
LEPTOGENESIS 

In this section we will show how an asymmetry is developed in the heavy neutrinos due 
to the curvature coupling term. This mechanism works for any model where lepton number 
violating interactions can take place freely at high temperatures. To be specific we will 
illustrate our mechanism using a specific GUT model like SO(IO). A recent example of 



leptogenesis due heavy neutrino decay with CP violation in a SO(IO) model is 2j|. In 
standard SO (10) unification, all Standard Model fermions of a given generation together 
with a right-handed neutrino are in a 16 representation of SO(IO), 

16f = (If + 5f + r0f)5t/(5) 

= (iV^ + (L,ci^) + (g,M^e^)) (25) 

The charged fermion and Dirac neutrino mass matrices receive contributions from Yukawa 
couplings of the form 16fl6fH (where H = IOh, 126h and/or 120h)- Majorana masses for 
the right-handed neutrinos are generated either from 



16f 16f 126h -DyS' N'i^ Nr (26) 

or from the non-renormalizable operators suppressed by some mass scale A 

j^WfWrlGHlGH^ J S^N^rNr. (27) 

When the GUT Higgs fields S' or S acquire a vev , a large Majorana mass M is generated 
for Nr which breaks lepton number spontaneously. This following the see-saw mechanism 
leads to small neutrino masses at low energies. At temperatures larger than the heavy 
neutrinos and the GUT Higgs masses one there will be helicity flip scattering interactions 
like S + Nr ^ S + N'^ which change the lepton number ( as T > M the helicity and the 



r{SNRi^SN'^) = {nsa) = ^(U T\ (28) 



chirality of A''^^ are same). The interaction rate is 

2 

The interactions decouple at a temperature To when TiTo) = H{To) from which we derive 
the decouphng temperature to be 

n^i3,.^[fj'^^^n,.^{^)'^^ (29) 

where we have used M = f{S)/A. Substituting fj29|) in the expression for lepton asymmetry 
( l24l ) we obtain the value of frozen in lepton asymmetry as 

Vl{To) = h7.8 X lO^^^^^s^^ (30) 

We choose the heavy neutrino mass M = lO^^GeV which is consistent with the atmospheric 
neutrino scale m^ = O.OSeV with a see-saw relation ttIj, = m'jj/M with the Dirac neutrino 
mass scale mo = 7GeV. Taking the GUT Higgs vev (S) = 1.2 x lO^^GeV and k^ = 0.07 we 
obtain the frozen-in asymmetry of the heavy neutrinos as 

7]l{Td) = 3.3 X 10-1° (31) 

We note that our scenario for generating heavy neutrino asymmetry is more general than 
the specific model we have considered above to calculate the temperature of decoupling of 
the N^ ^ N'^ interactions. We find that if this decoupling temperature is at the GUT scale, 

required heavy neutrino asymmetry 

-^-^^•^><^°( l.3xl0i^Gev ) (^^) (^2) 

This heavy neutrino asymmetry gets converted into light neutrino asymmetry when they 
decay, and light neutrino lepton asymmetry is converted into baryon asymmetry by the 
action of sphalerons in the electroweak era. 

The lepton asymmetry is passed on to the light neutrino sector when the heavy neutrino 
decays at temperature T r^ M r^ lO^'^^GeV. There exists AL = 2 interactions, vl + cjP <-)■ i'r + 



Td ~ 10^^ GeV then from (l24l) we see that if neutrino mass is M ~ 10^^ GeV we obtain the 



0° and VR+cf)'^* ^ z/j^+0°* that result from the effective operator Cw = '^{^La 0*)(0^ ^113) + 
h.c, where a, (3 denote the generation indices, /^o = {^a,Ga)J. ^^ ^^^ left-handed lepton 



doublet, (f) = (</>+, </>°)^ is the Higgs doublet {4> = ia2(l)* = {-(j)"^* , (p'f) [itJ. The The 
interaction rate for the interaction u^a + 0° ^ i^r/s + 0° is 

0.12 \Cap\'^T^ , . 

r = {n^ a) = — ^^^. (33) 

In the electroweak era, when the Higgs field in Cw acquires a vev, (0) = (0,^;)"^ (where 
V = 174 GeV), the five dimensional Weinberg operator gives rise to a neutrino mass matrix 
^™a/3 = "ivT^- Hence Cap can be written in terms of the heaviest neutrino mass rrii, as 
C ~ iTiyM/v'^. The lepton number violating interactions decouple when r(T;) = H{Ti) and 
the decoupling temperature Ti turns out to be 

Ti = la.Tvr V^ -^ ^ 2 10^^ f^_:^^) " GeV , (34) 

We see that the heavy neutrino decays occur at a temperature T c:^ M c:^ 10^^ GeV which 
is much below the temperature T/ = 2 x lO^^GeV at which the light-neutrino lepton number 
violating interactions are effective. Therefore the lepton number asymmetry from the decay 
of asymmetric number of heavy neutrino decays is not washed out by Higgs scattering with 
light neutrinos. 



5. CONCLUSIONS 

In conclusion, we have shown in this paper that a CP violating gravitational interac- 
tion between the heavy Majorana neutrinos and the Ricci curvature can generate a lepton 
asymmetry in thermal equilibrium in the GUT scale radiation era. The subsequent decays 
of these heavy neutrinos into the light standard model particles and the conversion of lep- 
ton asymmetry into baryon asymmetry can explain the observed baryon asymmetry of the 
universe. 

In order to generate lepton asymmetry we need two ingredients (a) a CP violating inter- 
acting term ([2]) of the curvature with heavy Majorana neutrinos and (b) a non-zero value 
of R in the radiation era to split the energy levels between N and N^. In [10| we used the 
warm inflation scenario to achieve a non-zero R at high temperature (we also studied the 
the role of beta functions of SU(N) gauge theories in generating a non-zero R). In this paper 
we have shown that a non-zero R can occur in the radiation era from loop corrections to 
the trace of the stress tensor and have calculated the contribution of the resulting curvature 

10 



perturbation to R. This mechanism of generating R is more general than the restricted 



(io|. 



scenarios considered earher in the earher study 

Some comments are however in order. In principle, one should also taken into account 
primordial perturbations of the gravitational background (characterized mainly by scalar and 
tensor perturbations) and of the energy density and pressure, characterized by 5p = 6Tq 



and 6p6l = 6T- (see for example |2l|). These perturbations are related as 6p = cj.6p, where 
Cg = w + pdw / dp is the adiabatic sound speed squared. For relativistic particles w = 1/3 and 
therefore c^ = 1/3. As a consequence, the trace of the perturbed energy- momentum tensor 
vanishes (this is not true in presence of anisotropic shear perturbations), so that according 
to the gravitational leptogenesis mechanism, no net baryon asymmetry can be generated. 

Moreover, we also compute the energy density of backreaction of quantum fields and 
compare it with the energy density of radiation. From Eqs. f lT2|) and flT3|) . and using 
H ~ l/2t, it follows 

(p) = {tI,^^'^) = 18A;i (2HH + H^ + WHH^) (35) 

, o; ij-4 3/C3 

The fci-term vanishes identically. Expressing the total energy density in terms of the scale 
factor a(t) = (aot)"*^'^ one gets 

P-.. + (P>-0 + ^. ^A,. (36) 

where pr = Tqq = -T^ is the energy density of the classical radiation. The ratio between 

oU 

the energy densities (p) and pr reads 



(P) fT 

r = 

where we have definite T* as 



4 



^-f-[¥j' P^) 



80 /15^^ 



3 



rCq 



For temperatures T < T* we have that r < 1, i.e. the energy density of quantum fields is 
subdominant with respect to the energy density of the radiation. This is in agreement with 
the approximation before discussed, where we pointed out that the effects of the backreaction 
wash out during the Universe expansion. In particular, since the decoupling temperature 
of heavy neutrinos occurs at GUT scales, T^ ~ lO^^GeV, we infer r ~ 10^^ ^ 1 and the 
backreaction is subdominant over the radiation density. 
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